ADA051896 


H 


mo  /fm  i 


AN  ALGORITHM  FOR  THE  ElZCTRjIlAG'KriC  SCATTERING  OUE  TO  AN 
AXIALLY  SYf^METRIC  BODY  WITH  AN  E-IPEDAiCE  BOUNDARY  CONDITini^  * a V 

A 


i 


r'ws'i 

! Ai’i-fCVvJ 

i DiaU 

1.  Dept,  of  Geophysics,  University  of  Utah 

2.  Dept,  of  Math.,  Old  Dominion  University 

3.  Dept,  of  Math,  University  of  British  Columbia  and  Dept,  of  Math., 
University  of  Utah 


•1  izi  pi.;blic  »l*ia3ej 
ihuduc.  Uniiniited 


Research  Supported  by  U.S.  Army  Research  Contract  #DAAG-29-76-G-0210 
at  the  University  of  Utah  and  by  NRC  Grants  A-6^01  and  A-8240  at  the 
University  of  British  Columbia. 


ii 

ABSTRACT 

Let  B be  a body  in  R^,  and  let  S denote  the  boundary  of  B. 

The  surface  S is  described  by  S - {(x.y.z)  : A + y ••  f(z),  -l^z^l  } 
where  f is  an  analytic  function  that  is  real  and  positive  on  (-1,1) 
and  f C±  1)  - 0.  An  algorithm  is  described  for  caiputing  the  scattered 
field  due  to  a plane  wave  incident  field,  under  Leontovich  boundary 
conditions.  The  Galerkin  method  of  solution  used  here  leads  to  a block 
diagonal  matrix  involving  2tl  + 1 blocks,  each  block  being  of  order  2(2N+1). 
If  e.g.  N - 00^),  the  conpated  scattered  field  is  aocurate  to  within 
an  error  bounded  by  Ce"*^**  where  C and  c are  positive  constants  depending 
only  on  f. 
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1.  IITTROD'JCriO:.'  AI'ID  S^I  u'^ARY 


Let  B be  abounded  body  in  R , having  surface  S vdvLch  is  given 


a.i) 


S “ {Oc,y,z)  : /x  +y  •*  fCz)»  > 


viiere  f is  cin  analytic  function  that  is  real  and  positive  on  (-1,1)  and 


Cl, 2)  a+z)^^  (1-z)^^  < fCz)  < a+z)®*  (1-z)®*.  -1<Z<1, 


where  , aj<l  and  6j<l  are  positive  constants.  In  this  paper  we 
describe  as  algorithm  for  ccraputing  the  field  scattered  from  B due  to 
an  incident  field  E° (r)  of  the  form 


a. 3) 


“ ik ‘r 
E (r)  » e e o 


where  e and  k (Ik  | = k =(i)/c  = lir/X)  denote  the  polarization  and 
0 0'  o 

^ ^ ^ ^ 

propogation  vectors  respectively,  and  r « xx  + yy  + zz,  where  x,  y 
and  z are  the  unit  vectors  pointing  in  the  direction  of  the  x,y  and  z 
axes  respectively. 

Let  the  body  B (rsp.  free  space)  be  homogeneous,  with  permit- 
tivity e (rsp.  e^)  permeability  y (rsp.  y^)  and  conductivity  a (rsp. 


a ),  so  that  the  refractive  index  of  the  body  is 
0 


where  w denotea  the  frequency  of  the  incident  field.  VJe  shall 
furthermore  assume  that 
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a. 5) 


iNl  Ik^l  p » 1. 


idiere  p is  the  smallest  radius  of  curvature  of  S.  This  assun^Jtion 
enables  us  to  apply  the  Leontovich  boundary  conditions  ( 15 , 21  ] 


a.6) 


(n  X E)  X n “ nZ  n X R 


on  the  surface  of  the  body,  where 


a.7) 


n - y/CPoN),  z “ (Vq 


and  where  n denotes  the  outward  unit  normal  to  S,  and  E and  H denote 
the  total  electric  and  magnetic  fields  on  S.  The  conditions  (1.5)  and 
(1.6)  are  satisfied  automatically  if  the  body  B is  perfectly  conducting, 
i.e.  if  a = <». 

The  condition  (1.6)  makes  it  possible  to  obtain  a singular 
vector  integral  equation  over  S for  the  surface  current  K on  S.  In  the 
present  {saper  vje  describe  an  algorithm  for  solving  this  integral 
equation  via  the  Galerkin  method,  using  as  basis  functions 

sinlTT  N*^{log(^^)-  ^ ) J » 

a.8)  (z,.)  £\z)  (1-zh 


m*  0,  ± 1,  ...»  - N,  n*  0,  ±1,  ±M. 


These  basis  functions  efie  jtively  VanJle  sinf^ularitios  of  K 
and  as  a function  of  z,  whidi  o=c’-r  at  z - +1;  they  are  siniilory  very 
effective  approximonts  in  the  (fr  variable,  since  the  Fourier  series  of 
^ (and  therefore  K).  converges  very  rapidly.  The  singularities  occuring 
in  the  kernel  of  the  integral  equation  for  K are  of  the  type  l/(z'  - z) 
or  log  I z'  - z]  at  z - z',  and  of  the  type  of  f°(z) , ra  ■ 0,  1,  ...  at 
z ■ ±1.  Ihe  first  of  these  is  effectively  handled  by  substracting  out 
the  principal  Vcilue.  The  renaining  ones  are  effectively  handled  by 
means  of  the  quadrature  formula  (see  [ 22  1 ) 


(1.9) 


r ^ 

g(t)dt  - h ^ 

2e'^^  g 

kh  , 
e -1 

J ^ k— L 

(1+e^)^ 

kh_^, 
[e  +1) 

after  transforming  the  intervals  (-l,z)  and  (z,l)  to  (-1,1). 

The  Integrations  over  S involve  two  variables,  4>'  and  z'  . While 

Che  integrations  with  respect  to  z'  must  be  carried  out  numerically,  due 
to  singularities  of  the  kernel  in  the  region  of  integration,  the  Integrations 
with  respect  so  (fi  are  carried  out  explicitly,  the  results  being  expressed 
via  hypergeometrlc  functions.  The  hypergeometrlc  functions  have  logarith- 
mic singularities  which  were  not  present  in  the  kernel;  for  this  reason 
explicit  integration  and  later  evaluation  of  the  hypergeometrlc  function  as 
described  in  the  Appendix  has  an  advantage  over  direct  numerical  integration, 
since  any  known  direct  numerical  integration  procedure  such  as  Che  trap- 
ezoidal rule  would  poorly  handle  this  type  of  transformation  of  singularity. 

The  use  of  the  basis  function  (1.8)  thus  leads  to  a block 

diagonal  Galerkin  system  of  equations,  one  syste.!!  of  order  2(2N+l)  for 
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each  m.  By  fo:Tning  2Mfl  such  blocks,  and  ta.king  N-M  , v.-e  arrive  an 
approjrimation  IC  of  K which  is  accui-ate  to  within  an  error  F,  inhere 
|el  ■ 0 (e"*^**  ) as  N -*■  » v;ith  c>0  and  independent  of  N.  The  use  of 
(1.8  ) furthermore  makes  it  possible  to  evaluate  the  scattered  field 
E*  by  means  of  simple  one-dimensional  trapezoidal  rule  integrations. 

The  error  in  our  approximation  of  tfie  scattered  field  E*  satisfies 


the  relation 


i 
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_ _ _ _ ^*5 

I'Sj^Cr)!  = |E®(r)  - E^(r)  1 = 0(e  ) as  N •*■  «, 

for  all  r on  the  exterior  of  B and  not  arbitrarily  close  to  B. 

The  above  problem  of  conputing  E®  given  E°  as  in  (1.2)  and  S 

as  in  Cl.l)  was  studied  in  [ 20  ] for  the  perfectly  conducting  case  and 

in  [ 5 1 for  the  case  as  described  above.  The  Galerkin  approximating 

basis  function  used  in  [ 5,20  ] are  of  the  form  -cos(ii»^)  S^(z)  and 

sln(iiip)  S (z),  where  S is  the  linear  spline  which  is  zero  at  z ^and 
n n ■ n-1 

2 

and  1 at  z^.  Thus  the  resulting  rate  of  convergence  is  0(1/N  ) 
if  f has  no  singularities  at  z * ±1  and  0 (lyil^)  if  e.g.  f(z)  ~ Cd-z)*^ 
as  z-^1,  where  it  is  assumed  that  the  interval  (-1,1)  is  divided  into 
N subintervals.  In  addition,  the  qxHdratures  used  in  I 5,20  J converge 
very  slotvly  as  a result  of  the  singularities  present  in  the  integral 
equation.  Furthennore  the  expression  for  the  gradient  of 
0 . elH.  II-' 'I /(%  I r-r' |)  obtained  in  [ 5,20  ] is  incorrect. 

The  algorithm  of  the  present  paper  has  been  checked  out  on 
a computer  for  the  case  of  a sphere  of  radius  1 for  which  the  surface 
current  K and  the  scattered  field  ^ can  be  expressed  explicitly. 

Using  M - 4 and  N=M  = 6 , we  find  that  K is  accurate  to  4 significant 
figures,  and  for  r > 2,  E®  is  also  accurate  to  4 significant  figures. 
The  paper  is  organized  as  follows. 

In  Sec.  2 we  describe  the  geometry  of  the  surface.  In  Sec.  3 
we  derive  a representation  on  the  surface  S for  the  incident  plane  wave. 
Sec.  4 contains  a derivation  of  the  integral  equation  for  the  surface 
current,  as  well  as  an  integral  expression  for  the  scattered  field  in 
terms  of  this  surface  current.  In  Sec.  5 we  describe  the  basis  functions 


a. 10) 
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to  be  used  in  the  Galerkin  r^ethoc  cf  Sec.  7.  Sec.  6 contains  an 
approximate  representation  of  the  incident  electric  field,  in  terms 
of  the  basis  functicxis  of  Sec.  5.  In  Sec.  7 derive  the  Galerkin 
equations  for  the  surface  current,  and  we  describe  a method  of  computing 
the  coeffieicents  of  this  system,  and  for  solving  this  system.  In  Sec.  8 
we  describe  a procedure  for  evaluating  the  scattered  field.  In  Sec.  9 
we  discuss  the  rate  convergence  of  the  procedure.  Finally,  Sec.  10 
contains  a numerical  example,  illustrating  the  algorithm.  Appendix  A 
contains  a study  of  the  functions  derived  in  Sec.  7 as  well  as  of 
their  derivatives.  The  results  of  this  appendix  illustrate  the  type 
of  singular  behavior  of  the  functions  and  thus  they  dictate  the 
type  of  approximate  methods  to  be  used  in  order  to  achieve  high  accuracy, 
and  they  sin^jlify  our  proof  of  convergence. 

The  rate  of  convergence  of  the  method  of  this  paper,  namely 
0(e  ) using  an  approximation  of  the  form 


a. 11) 


M 

I 


N 

I 


m=-^  n=-N 


a 0 (z)e 
mn  n 


lia.9 


for  each  ccmponent  of  the  surface  current  is  optimal,  in  a certain  sense. 
By  the  results  of  [ 25  ] , given  any  approximation  method  of  the  type 
Cl. 11)  which  is  to  converge  for  all  f analytic  on  (-1.1)  and  satisfying 
Cl.2),  the  resulting  error  of  this  method  cannot  converge  to  zero 

u 

-yN  * 

faster  than  e ' , for  some  y > 0. 

Similar  numerical  methods  have  been  considered  but  it  is  believed 
that  the  order  of  convergence  obtained  is  not  as  good  as  that  demon- 
strated here.  For  further  discussion  of  such  methods  see  for  exanple 
Andreasen  [ 2 1 whose  proposed  method  considers  a maxinum  period  of  20 
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wave  lengths.  Barber  and  Yeh  [ 3 1 and  V.'aterrran  (27  1 have  considered 
extended  boundary  jnsthods,  while  Kennaugh  [ 13  1 and  Schultz,  et.  al. 

[ 19  1 have  discussed  other  inplementations  using  a product  z,  <))  basis. 
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2.  GEOMETRY  OF  THE  S'JRFACE 


A point  on  the  surface  S is  represented  by 


C2.1) 


r = f(z)  cos<fi  'X  + fCz)  sin^  y + zz. 


Ak  A 


vAiere  x,  y and  z denote  the  unit  vectors  in  the  direction  of  the  x,  y 
and  z axes  respectively. 

It  is  convenient  to  introduce  three  unit  vectors  on  the  surface. 


A A A 


n,  ifi  and  t where; 


C2.2) 


n a a(z)  cos  ifix  + a(z)  sinv>  y - f Cz)a(z)z 

AAA 

=~sin'^  X + cosV  y 

A A A A 

t a f'(z)o(z)  cosv^  X + f'(z)a(z)  sin^  y + a(z)z. 


where; 


C2.3) 


a(z)  a 1 1 + f"  (z)^  1 ^ 


The  vector  n is  the  unit  nonral  to  the  surface,  0 is  the  unit  vector 

A 

at  r,  pointing  in  the  direction  of  increasing  and  t is  the  unit 
longitudinal  vector,  pointing  in  the  direction  of  increasing  arc  length. 

AAA  AAA 

Thus  n,  ip,  t and  x,  y,  z are  related  by  means  of  the  equations 


r 


j 

S'! 


f: 

i 


(2.4) 


and 


(2.5) 


V 

/ 

n 

v> 

“ 

tj 

1 

0(2)  cos.^ 
-sin'/> 


a(2)sin‘^ 

cosi^ 

cosV’  £'(2)0(2) 


-f  (2)0(2) 
0 

sinV>  0(2) 


X 

y 

A 

'a' 

X 

o(z)cosip 

-sinip 

A 

y 

s 

o(z)sin'P 

COS'i’ 

A 

z. 

-f  (2)o(z) 

0 

f (2)o(z)cos';> 
f ' (z)o(z)siniP 
a(z) 


n 

A 

A 

J 

.t 

I 

i 


* ^ 


C3.1) 


3.  TIE  INCIDEIT  ELECTRIC  HELD 


Let  the  incident  radiation  be  a plane  vave,  given  by 


r.o,— , — ikrt'r 

E Cr)  = e e o , 


vdiere  e points  in  the  direction  of  polarization,  and  is  the  pro- 
pogation  vector  which  satisfies  the  relations 


C3.2) 


, - ir  1 2ir 

k = k = — = -r- 
O O C A 


k • e = 0 
o 


We  shall  furthemvore  assume  that  k lies  in  the  xz  plane 

o 

and  makes  an  angle  0^  with  the  z -axis.  Thus 


(3.3) 


k = k sin  0 x + k cos  z. 
o o 0 o 0 


It  is  convenient  to  set 


(3.4) 


where  and  a2  are  scalars,  vdiile 


(3.5) 


e,  “ y,  = -X  cos  0 + z sin  0 

I » ' y o t 


1 
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Borison  [ 4 ] has  shovm  that  if  a2  = 0 (rsp.  = 0)  then  the  back 
scattered  electric  field  E®Cr)  is  polarized  only  in  the  direction 
(rsp.  C2). 

Using  (2.5)  and  (3.5)  we  can  express  F in  components  of  s?,  t 

A 

and  n.  We  get 


e =»  a^[  a(z)  n + cossp  s?  + a(2)  £'(2)  sinv’  t ) 

(3.6)  + aJ-<x(z)(cos  0 coav*  + sin  9 f'(z))  n + cos  0 slxvfi  <p  ' 

Z 0 0 0 

-a(z)(cos  0^  cos  f’(z)  - sin  0^)  t ] . 

Next,  let  us  find  the  Fourier  expansion  of  E°(r)  on  S.  To  this 
end  we  use  the  identity 


/-  r7s  Ixcosv’  r .Jn  . . Imp 

(3.7)  e = > i J (x)e 


where  J (x)  denotes  the  Bessel  iiuiction,  J (x) 
m m 

Using  (2.1)  and  (3.3)  we  get 


” (-1)^  (x^/4)°^^ 
(n-hn).' 


(3.8)  k • r = k f(z)  sin  0^  cosip  + k z cos  0 • 

O O 0 O 0 

Hence  combining  (3.1),  (3.7)  and  (3.8)  we  find  that  the  incident  field 

on  S is  given  by 

(3.9)  E°(7)  - F *0  ? l"j  (fc  fCz)sln8  )e‘™’. 

-moo 

r3=-" 

where  e is  given  in  (3.6).  Combining  (3.6)  and  (3.9)  we  get 


iM 
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^ ''1 


A 

E°(r)  « {la^a(z)  sinv>  - a^cxCz)  (cos  6^  cosv>  + f’(z)  sinS^)  In 

A 


(3.10) 


+ I a^cos<P  + &2  ®o  ^ 

+ I aj^a(z)f ' (z)  - a20((z)  (cos  0^  coss?f'(z)  - sin  0^)  1 t} 
•.‘"o'  ®0  1 i”j„(k„t(0  Bio 

QJS— 00 


[ S 
t f 


H 


The  scattered  field  E®  - E®(r  ) is  given  in  terras  of  the  total 
electric  (E)  and  magnetic  (H)  fields  on  S by  [26] 


(4.1) 


E® 


I iu)y(nxH)G  + (nxE)x7G  + (n*E)VG  ]dS 


where  it  is  assumed  that  the  field  vectors  E and  H have  time  dependence 
of  the  factored  form  e and 


(4.2) 


G(r,r')  - exp  [llc^|r-r'|  ] 


r-r' 


The  renaining  vectors  in  the  integrand  (4,1)  are  expressed  in  terms  of 
7,  and  V is  expressed  in  terms  of  7. 

Let  K denote  the  surface  current  on  S,  due  to  the  fields  E and 
H.  In  view  of  the  Leontovich  boundary  condition  (1.4),  K satisfies 
the  relations 


(4.3) 


K 


- nxH  ; nxE  - - nZ  nxK 


a»  - en-E  > ^ V-K 


Our  development  of  the  integral  equation  for  K follows  that  in 


I 5 1 . We  include  the  derivation  in  ( 5 ] for  this  equation,  for  sake  of 
ooftpleteness. 
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Let  A be  a continuous  vector  field  tangent  to  S.  Then  the 
following  results  are  valid,  if  S satisfies  the  Lyapunov  conditions  ( 17, 
p.  90  1 (a)  the  integral 


I(r')  5 A(r)  G(r,r')  dS 


Is  a continuous  function  of  r'  In  R"’. 


(b)  As  r'-*T'  4 S , the  relations 
o 


n(r  ) X lira  A(r)  x VG  (r,r’)dS 

° r'-*T'  - 

r o S 


±H  A(r')  + n(r  ) x lA(r)  x VG  (r^,r)  JdS 
o o o 

Js 


are  satisfied  where  the  plus  (resp.  minus)  sign  corresponds  to  an  approach 
from  the  outside  (resp.  inside)  of  B. 

(c)  The  term  [17,  p.  95  ] 


n(r')  • A(r)  x VG  (r,r')dS 


is  a continuous  function  of  r'  on  S.  The  term 


E^Cr’)  H n(r)  • E(r)  VG  (r,r')dS 


suffers  a discontinuity  on  transition  through  S equal  to  n*AEj, 
vrtiere  is  the  difference  of  the  values  outside  and  inside . The  third 
term  in  E®,  therefore  does  not  affect  the  tangential  component , but 
reduces  the  normal  conponent  of  E to  zero. 

Since  the  total  electric  field  E - E°  + E*  is  zero  inside  the 
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scattervar,  (4.1)  yields 


0 « n ^ i°  - _ Um  n(r')  ( iu)y(nxH)G 

r'-^'eS  ° 

° Js 


+ (n?«E)  X VG  + (n*E)  VG  1 dS 


where  the  approach  is  from  the  inside. 
Application  of  (b)  gives 


C4.4) 


0 - n(r')  X E°  (r')  + n (r')xE(r') 

- n(r')  X f I la)p(nxH)G  + (nxE)  x VG  + (n*E)VG  ) dS 


Next,  tcdcing  the  limit  as  jr*’r'  from  the  outside  of  B in  (4.1) 


and  using  the  relation 


n(r')  X E(r')  " x ■*"  (r')) 


we  get 


nxE  = nxE°  - _lim  n(r’)  x I lwy(nxH)G  ] 


r'-^r' 


+ (nxE)  X VG  + (n*E)  VG  ] dS 


Application  of  (b)  yields 


(4.5)  I n(r’)  X E(r')  - n(r')  X E°(r)  - n (r')  X I iup(nxH)G  + 

+ (nxE)  X VG  + (n*E)VC  1 dS 


T?1E  a'^SI  S FUNCTION'S  FO;-<  APPROXIMATIN 
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Let  u;3  set 


(5.3) 


v(z)  - f(z)  (1-z^) 


vd>ere  f e satisfies  (1.2),  f # 0 in  (),.  Let  H - H (d,d')  denote 

a u 

the  fandly  of  all  functions  F ■ F(z,w)  such  that 

i)  F(z,e^''’)  / v(z)  belongs  to  H (fl^)  as  a function  of  z for 
all  v>e  I 0,2ff  ] ; 

ii)  F(z,w)  belongs  to  H(Aj,)as  a function  of  w for  all  ze  (-1,1  ] . 


We  define  a norm  on  H(d,d')  by 


(S.4) 


If  F e H(d,d'),  we  approocinate  F on  S ■ (-1,1  1 x i C,2it  ] as 


follows 


(5.5) 


S (2,»)  5 ! <*•*’) 

n*»-N 


where  the  a are  numbers  and  the  (j)  are  basis  functions  given  by 

ma  Titti 


♦i»  <'•'>  = 

0 (z)  = 1%)  d-z')  sine  I 


(5.6) 


wCz)  = log  (^I^)  ; sine 
h > 0. 


8ln(irx) 
* “ TTX 


The  numbers  a are  given  by 

nn 


18 


(5.7) 


a 

Dm 


2Trf^(z  )(l-z  h 
n n 


2ti 

ill?. 

F(z  ,e  )e 
n 


z ■ tanh  (nh) 
n 2 


d'^ 


Next,  recalling  the  definition  of  f and  the  relation  (1.2),  let  us 


set 


(5.8) 


Y - >5min  (a2»S2^  ’ ^ (irdY2)?log  d'  ] 


).l:  Let  F s H,  let  N •»  M , and  let  h - [ ird/  (Y2N)  1 • Then 
lliere  exist  constants  C,  and  C2  which  are  independent  of  N,  such  that 


(5.9) 


I (».*’)  I < CdV^'"' 


(5.10) 


”*«  ' i F(z.«)  - i (z.F)  1 < Cj»e 


(z.v»)es  I 3z 


(5.11)  I i - I V I ^ "2"  ^ 


-yn" 

h 


H -YN 


Proof:  It  is  shown  in  [ 16  ] that  if  g/v  g H (O^j)*  then  by  taking 
h “ (■ird/(Y2N)  there  are  positive  constants  C'  and  C"  such  that 


(5 


N f^Cz  ) _yN 

12)  I r(z)  - y “7“?  9 (z)  < C’e  ^ 


and 


(5.13)  , I *’(z)  - ! 4^9;  (»)  I 

zQ-1,11  n— N " 


19 


Similarly,  it  follows  from  Cauchy's  theorem,  that  if  H(A.,), 


and  if  a is  defined  by 
m ^ 


(5.14) 


m 2Tr 


2it 


G(e^'^)  dv»,  m— M,-Mfl, . . ,M, 


then  there  are  constants  K'  and  K"  such  that  for  all  ^6  I 0,2ir  ] , 


(5.15) 


M 


I GCe^'^)  - I a 1 < K'  (d')"“ 

__  X.  m “ 


m— M 


and 


(5.16) 


M 


||.  I iBay*  I <K"  (d’)-" 

np-M 


m 


On  noting  that  M - N , the  inequality  (5.9)  is  obtained  if  we 
use  (5.11)  and  (5.15)  in  thm  6.2  in  ( 24  ] . Similarly,  (5.10)  is  a 
consequence  of  (5.13)  and  (5.15),  while  (5.11)  follows  frcm  (5.12)  and 
(5.16). 


>2. 


Theoren  5.2  [16  ] . If  vg/f  *6  ®(^)^  » then  there  exists  a constant  K such  that 


(5.17) 


-1 


V(z  )g(z  ) 

g(z)  9 (z)  dz  - h \ ^ I < j,  -TTd/h 

f *(z  .)a)'(z  ) 
n n 


Tliis  theorem  shows  that  if  g is  any  function  such  that 
vg/f^H(ftj),  then  for  h sufficienty  snail,  the  sequence  {0  } may  be 
considered  to  be  an  orthogonal  sequence,  for  practical  purposes,  and 
g may  be  expanded  with  respect  to  this  sequence.  The  coefficients  of 
this  expansdlon  take  on  the  very  sriple  form 


■ : 
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I. 


(5.18) 


f *(z  ) 


tanh  (nh/2) , 


For  purposes  of  numerical  integration,  we  state  the  following. 


on 

i. 


Theorem  5.3  I 23  ] : Let  g 6 H , and  let  |g(z)  | < C (l-z^)“"^ 
(-1,1),  where  a > 0.  If  N > 0 is  an  integer,  and  h = (2ird/aN)’^, 

then  there  exists  a constant  C’ , independent  of  N,  such  that 
fl 


(5.19) 


N 8(z_) 

sW^-h  I HTTf-) 

n=-N  n 

-1 


< C e 


-(2TrdaN) 


Finally,  thra.  5.4  which  follows  describes  the  accuracy  of  the 
midordinate  rule  used  in  Secs.  7 and  8 for  integrating  periodic  functions. 
Theorem  5.4:  Let  d'  > 0 and  let  kG  H(A,,)  be  bounded  on  A.,.  Then 

— — <1  d 

there  exists  a constant  C depending  only  on  g such  that  for  M = 1,2,3,..., 
2Tr 

1^, 


(5.20) 


St*"')  I 

k=»l 
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6.  APPROXIMATION  OF  TtlE  INCIDElNiT  ELECTRIC  FIELD 
We  shall  make  the  approximation 

_ _ M N ^ _ 

(6.1)  E°(r)  : E°(r)  = J I I a ' ^ + 6’  t ] ^l>.r  e s, 

**  ®n  nm  inn 
m»-M  n*-N 

of  the  incident  electric  field,  where  il;  are  defined  in  Sec.  5.  For 

’ ^mn 

the  sake  of  convenience  we  shall  use  the  notations 


z “ tanh  (nh/2) 
n 

(6.2) 

“ 1°  exp  (Ik^z^  cos  9^}j^(k^f(z^)8lii  e^) 


As  will  be  seen  in  Sec.  7,  we  will  approximate  a sligjitly  altered  field 
j®  = vE°  where  v(z)  ■ f(z)  (1-z^).  Tnus 

_ M N 

j’’(r)  “ J°  = ^ ^ [ a + 6 t 1 . 

^„  inn  im  inn 

in»-M  n*-N 

The  results  of  Secs.  3 and  5 show  that  if  r e s then  each 
coirponent  of  is  in  H(d,d').  Using  the  formulas  (5.6),  we  have 


(6.3) 


. f*(z  ) 

a = . n 

"IF' 


6 - 

mn  ' 


2tt 


2tt 

F’  (r)  • s?  e ds?  | 


0 

2ti 


E^(r)  • t e d>fi 


z*z 


z*zv 
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7. (a)  THE  GALERKIN  EQUATION  FOR  THZ  SURFACE  CURRnTT 
DERIVATION  OF  THE  EQU-AriONS 


Rather  than  solve  (4.11)  for  K,  it  is  numerically  more  convenient 


to  define  J by 


(7.1) 


where 


J(r)  . v(2)  K(r) 


(7.2) 


V(2)  - f (z)  (1-2^), 


and  to  solve  the  resulting  integral  equation  for  J.  This  transformation 
helps  to  take  account  of  the  unknown”  singular  behavior  of  K at  z = ±1, 
and  enables  us  to  effectively  approximate  both  J and  its  first  derivative 
with  respect  to  z,  by  methods  of  Sec.  5. 

The  substitution  (7.1)  replaces  (4.  6 ) by  the  equation 


(7.2) 


{isnzJ+v  J G + (rrcJ)  XVG 


+ ~ (V-ij)  VG  ] dS  + J°  } = 0 

tan 


vdiere 


(7.3) 


*IF  S satisfies  Liapunov  conditions  (see  Sec.  4),  then  K is  bounded  on  S. 

It  is  difficult  to  determine  the  exact  singular  behavior  of  the  derivatives 
of  K at  2 ■ ±1. 


i 
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We  now  make  the  Galerkiia  ajor'cxl’nation 


(7.4) 


_ M N 

J ■ J “ I I I a I?  + b t ] -<1) 

**  ran  ran  ran 


m=-M  n=-N 


in  (7.2),  vdiere  the  \h  are  defined  as  in  (5.16),  and  the  a and  b 

’ ’^ran  ran  ran 

are  unknown  numbers. 

Let  us  now  recall  that  if  7,  r ' e S,  then 


(7.5) 


R = |r-r'|  = - 2ff*  cos  + (z-z')^} 

, ik  R 

r = i e ° 

® 4tt  R 


where  here  and  henceforth 


(7.6) 


f = f(z)  , f*  = f(z')  , a ■=  a(z),  a*  - a(z'), 


We  shall  also  use  the  notations 


(7.7) 


3z 


f ^ f-f*cos(sg-»t> ') 
“ 3f  “ R 

z 3R 


z-z 


ifl  ff*  sin 
^ “ R 

(G^.G^.C'’)  - (R^,r'',r'^)  ^ G. 


M N 

Moreover,  v/riting  I for  I ^ , 

m,  n m“-M  n“-N 


vje  have 
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Substituting  (7.4),  (7.10),  (7.11)  as  v:ell  as  the  approximation 


(7.12) 


S = I I a v»  + B t ] 

^ Oin  DD 


into  (7.2),  and  recalling  that 


(7.13) 


dS  - i d;>  dz. 


we  get 


(7.14) 


{ %!  Z J [ a_v>*  + t*  1 ij)* 

**  mn  im  nni 

in,n 

rl  fZlT 

+ V*  Gy(a^  + b t]\j) 

ov  ^ nn  mn  ^nm 
m,n 

Jo 

+ ^ J I a f(G^-f'G^)i  + b^  f(G^-f’G*)  t 

\j  Tfin  mn 


V mn 

m,n 


- { a - G*^  + b f(f'G^+G*)  } n ] 

mn  a mn  mn 


+ { i-  g'^  i + (f ’G^+G^)  t + (G^-f ’G^)  n } 

lue  oc 


y [ fsHL  + a b l-r  C-  'I'  ) 1 } dz 

i ‘ V mn  3z'  S)  ^mn' 

m,n 


+ y ( a <fi*  + & t*  1 <>*  } “ 0 

^ * rnn  mn  inn 


where  the  starred  variables  denote  functions  of  z'  and 


I 
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It  is  convenient  to  introduce  several  identities  in  order  to 
reduce  C7.14)  to  a system  of  linear  algebrciic  equations.  The 
equations  (2.4)  yield  the  identities 


(7.15) 


• t* 


cos  (v>-^'  ) 

- a(z)  f'(z)  sin  (s?'-v) 

- a(z)  sin 

a(z')  f'(z')  sin  (v>'-v) 


t • t*  ■ a(z')  a(z)  I 1 + f ' (z)  f'(z')  cos  (^-v')  ] 
n • t*  - a(z')  a(z)  [ f'(z')  cos  - f'(z)  ] , 


These  identities  are  useful  for  takii'g  oonponents  of  <fi*  and  t*  in  (7.14). 

The  identities  (7.16)  - (7.23)  which  follow  serve  to  achieve 
further  sinplicity  by  enabling  us  to  symbolically  elimii'vate  the 
integrations  with  respect  to  . 

Upon  setting 


(7.16) 


ik  r-r'  " 
e o 


G(r,r')  - = I 

4ir  r-r' 


im(v>'-V>) 


it  follows  that 


(7.17) 


.i'*’’  c - 

m 2t7 


2TT 

G(r,r')  di^  . 


Notice  furthermore,  that  - G^.  Therefore 


J 


1 iiiV>' 

(7.18)  27r  G(r,r')  cos(^^'-V>)  dss  - 2 ^ ^ ^n-1  ^ 


(7.19)  ^ G(r,r')  siu  («^-v')  ^ 2i  ^ ^ntf  1 ” ®m-l  ^ 


By  means  of  integration  by  parts,  we  furthennore  find  that 


(7.20) 


(7.21) 


(7.22) 


1 9G(r.r').  , ziS  G(r,r') 

W ^ 2 TT 


. iniP'  „ 

- im  e G„, 
m 


ijj-  e^™’  cos  (v>-^')  |^dv> 


i [ (iiH-1)  + (m-l)  1 > 


1 ira^  . 9G(r ,r * , 

^ e sin  (-^-^  ) 3^ 


I ("H-D  G^,  - («-l)  G 1 . 


We  use  C7.13)  in  (7.12)  to  take  conponsnts  of  and  t.  Then, 


recalling  that 


(7.23) 


(z,^)  - e (z) 

mn  n 


vdiere 


C7.24) 


e M - sine  I 1 

" f(s)^  '• 


and  where  vfz)  is  given  in  (7.2),  and 


(7.25) 


- log  (^). 


we  can  synibolically  carry  out  the  integrations  with  respect  to  in 
the  equation  resulting  from  (7.12),  by  using  (7.16)  - (7.21).  Upon 
equating  coefficients  of  in  the  resulting  equations,  we  obtain 


(7.26) 


Jsnz  [a  0 +vItP™a_  + Q"”b^l  - 
^mnn  mn  mn 

n n--N 

Jsnzyb  0 + vy  [ R™  a + b ] 

^ ^ mn  n ^ „ nin  inn 

n n— N 


“ M,  “Mtl  ,...,  M. 


- y o e 

^ mn  n 

n^.ff 

N 

- - y B 0 

*'  mn  n 

n«-rt 


The  coefficients  p"“,  q”",  r'^  a.nd  s““  depend  on  z',  and  are 


given  by 


:7.29) 


(7.30) 


Next, 


..on  _ 
R -TV 


_n  ( G _ I + i(itrf-l)n  Z+  1 


nri-1 


-1 


+ G , a*f* 
m-l 


M ..  z - 1 

+ G a*  I - 2inf'  n Z 1 

IQ 

+ G^  a*f*'  [ 1 


! 


+ dz 


u)ea 


o“» 

S “IT 


E • -J  /.«!  f f * 


(^){  a*f*'  W ff  I + Vl  ^ 


-1 


i 


^ ('>  < ‘ Vi  - Vi  * 


lue  m 


dz 


I Z 

In  these  equations  f*'  - f ' (z' ) . «*  “ I ^ ^ ^ 


«ttiJig  z'=  - tanh  Uh/Z)  in  (7.26)  and  using  the  relations 


0 if 


v(z^) 


(7.31) 


if  n-  I 
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we  arrive  at  the  system 


N 


Z a . + Cz  ) I [ P a + Q b 1 - - a 

^ A ^ * « mn  ^ 6 tnn 


(7.32) 


n»-N 
N 


£ nn  £ mn 


iq£ 


■i"  2 + f**  (..)  I I - S, 


n»-N 


nn  £ nsi 


in£ 


where  we  have  used  the  notation 


pim  _ pina 


C7.33) 


m „inn  , . _mn  „nn  , . 


The  systan  (7.32)  is  a block  diagonal  system  of  the  fom 


(7.34) 


-M 


-Mfl 


'M 


1 

^ 

1 — 

I 

• 

• 

• 

m 

“->H-l 

• 

• 

• 

1 

P 

3C 

1 

where  each  is  a conplex  matrix  of  order  2(2N+1),  ard  (since  G_^  “ G^) 

B ■ B . The  m'th  system 
m -m  ■' 


(7.35) 


B a “ a 
mm  m 


in  (7.33)  corresponds  to  all  of  the  equations  (7.35),  for  fined  m. 

TRU  • 

Thus  if  we  denote  by  A the  2x2  matrix 

X 


34 


7.(b)  EVAUJATION  OF  P™.  R°"  and  S°” 

It  is  ccxivenient  to  set 


(7.39) 


8 /V 
n 


0 f 

(_!!_)  t 
' V ' 


6 

n 


The  following  integrals  app>ear  in  (7.27)  - (7.30). 
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1 


it/a)  G dz 

IB 


Y ff  G dz 
n m 


1 


G dz 
m 


1 

Y f G 
n m 


I Y„/(af)  1 G„  dz  r 


Y f G dz 
n m 


Y f G"  dz 
'n  m 


(7.40)  J' 


Y„  (f/o)  G"  dz 
n D 


6 (1/f)  G dz  (lai’K)) 

n m 


6 G^  dz 
n m 


6 G‘  dz 
n m 


Y ff’  G dz 
n m 


(1/a)  G^  dz  (hp^O)  - Y„  (l/a)  G^  dz 


(f/«)  G^  dz 


Die  notations  (7.40)  enable  us  to  express  (7.29)  - (7.30)  in 
the  "mcare  conputable"  form 


P»  . , i i.M  - nz 


<»-«I  "In  I 


(7.41) 


'“^nrt-l.n  '^m-l,n  iirH,n  m-l,n 


+ _2_  (j(2)  _ j(2)  ) } 


(7.42) 


q“  . n { ™ - f-  CK^>_„  - K<«  „) 

- s?  1 '•^l.n  ■"  ■"  “^1.”  ■ ’ 


, „.£»■  {-  - l“>  „)  + inz  1 (irt-l)  + <"-l)  I,«.l,„ 


(7.43) 


n,  . (3)  1 21nZm  ,(5) 

+ Vl.n  " Vl.n  ^ ■ f*’  mn 

22.  ) 

i '^‘^uri-l.n  m-l.n  nrH.n  m-l,n 

+ _HL(j(2)  + j(2)  + 2 k^))} 

(ije  ' nri-l,n  in“l,n  f*  nm 


rr  ■ ^ 


} ' 

ft  : 
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y 

.v' 

.1 

-( 


(7.44) 


s“"  = ir  a*f*'  { iuy  (1^,'  + I 


(4)  ^ ^(4) 


T li«) 


nfl.n  m-l.n  f*  mn 

■ ^ n n " ^ 

iD+l.n  m-l,n  r*  mn 

= + I (m+l)  l1!J  ^ - (m-1)  + 

Itoe  nr(-l,n  m-l,n 


>, 

j-' 


nri-l,n  m-l,n  f*  mn 


V. 

i, 


The  integrals  (7.40)  can  be  sinultaneously  evaluated  for  all 

m,n  by  evaluating  the  quantities  f,  f,  a = (l+f^)  f*'  and  a*  = (1+f*'^) 

0 (n  - - N,  - N+1,  ....  N)  .and  G (m  » - M,  -m N).  The  integrals 

n m 

and  are  evaluated  sinultaneously  by  means  of  the  formulas 
mn  mn 


(7.45) 


tl 

•z ' 

H (z,z’ ) dz  “ 

H (z,z')  dz  + 

4 

ri 


H (z,z')  dz 


-1 


-I 


^ ^ u H (z’.y^)  + I H (z',Vj.) 


The  definitions  for  the  exact  form  of  u and  v follow  in  (7,48). 

8 S 

On  the  other  hand  the  integrals  K^,  and  are  evaluated  by 

means  of  the  following  fomulas.  i 

j , II  Here  vje  use  the  results  (A.2). 
mn  mn 


(z,z’)  g (z)  dz 


m 


-1 


(7.46) 


j - - -ij  • f*a*^  g(z')  log  (^:;|r)  [t 


4ir 


(G^  (z.z')  g(z)  + * |r"  z-z*  ^ 


-1 


r' 

> i 


1 ? 
c 


K , N . Using  CA.2), 
mn  mn 


1 2 
gI  Cz.z')  g(z)  dz  “ - 8(z')  log  (^) 


m 


4t 


(7.47)  -1 


1 2 

[G^  (z.z')  g(z)  J 

“ Air 


-1 

and  the  integral  on  the  right  hand  side  is  again  evaluated  by  means  of 
(7.45). 

We  split  the  integration  in  (7.47)  since  (ztz')  has  a 
singularity  at  z = z'.  The  formula  (5.18)  is  then  used  to  evaluate 
each  integral.  Thus 


(7.48) 


y'+l 

u ••  — h*  sech 
s 2 


l-z*  , . ,2  i 

= — 2~  h*  sech  > , 


1+zL  ^ tanh  (^) 


+ 2 


,th?. 

2 


for  suitably  chosen  h*  > 0. 


8.  APPROXU-IATION  OF  THE  SCATTERED  FIELD 


The  scattered  field  E®  is  expressed  in  terms  of  K by  means  of 
the  integi’al  (4.  7 )•  Once  T has  been  obtained  by  means  of  Secs.  7,8 
and  9 , we  form  an  approximation  K of  K by  means  of  the  equation  (see 
Eq.  (7.1)) 

(8.1)  K » K (z.^P)  = J (z,¥>)» 

and  we  substitute  K for  K in  (4.12)  to  get  an  expression  for  an  approxi- 
mation E®  of  E® . In  this  section  we  shall  give  a detailed  description 
of  the  evaluation  of  E®. 

We  shall  approximate  E^(r)  for 

(8.2)  7 ' = p cossp'  + p sins?'  y*  + z'z’ 

where  p > f(z').  If  p is  close  to  f(z)  i.e.,  if  7'  is  close  to  the  sur- 
face, (say  l7'-s|  ^.2  if  N = 10)  then  we  recoinnend  the  integration 
methods  of  Sec.  9 to  evaluate  the  integrals.  This  would  involve 
splitting  the  integrals  from  -1  to  1 into  integrals  from  -1  to  z'  and 
from  z'  to  1,  as  in  (7.45).  For  sake  of  sijiplicity,  we  shall  describe 
an  algorithm  for  evaluating  E*  which  is  valid  if  r'  is  not  arbitrarily 
close  to  S (say  |7'-s|  > .2  if  N > 10).  In  this  latter  case  it  is 
convenient  to  integrate  by  parts  in  the  integrals  with  respect  to  z 
which  involve  0'  , so  that  the  resulting  integrals  involve  6 . This 


» r.  -jm 


r ^ 


(8.6) 
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.3  ~ 


I I 

m,n 


_Rn  . . -inn  , 

a R + b S J 

xim  mn 


\ 
\ 

. i 


. *?■ 

i V 


(8.7) 


,s  ~ 


y I a t“"  + b U®*  1 
^ nm  nm 

m,n 


The  relations  (7.7)  - (7.11)  and  (7.15)  - (7.22)  enable  us  to 
obtciin  e:qplicit  expression  for  P®',  ...  , U®'.  Setting 


(8.8) 


9n(2) 

“n  " ‘^n  ’ 


we  get 
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n2iT  . 

(if  fcl, 

-10 


I a (-  sinv;  X + cos^  y) 
mn 


+ b (fa  cosv>  X + f'  a sini^  y + a z)  ] e o) 
inn  Q 

- ^ ^ 

+ nZ  I la  f(G*-f'G^)  (-  sinv>  x + cos^  y) 


(8.9) 


^ A ^ A 

+ b f (G  -f 'G^)  (f*  a cos^  X + f ' a sin^  y + az) 
nm 

- { a — g"^  (a  cosifi  x + a sin<^  y - faz) 
uoi  a 


+ b f(f'G^  + G^)  (a  cos*^  X + a sln^  y - f'az)}  ] e'^“'^w 

Trtn  u 


A ^ 


- /S  ^ r\ 

+ (f'G*  + G^)  (f ' a coss?  X + f a elw  y + a z ) 

+ (g‘  - f’G^)  (a  coss»  X + a sin(P  y - faz)  } 


* r (a  im  0)  + a b (foJ  ) I e } dv>  dz 
*-  inn  n nm  n 

m,n 

i.k  R* 

Here  G - ; R » { (z-z')^  + - 2f  pens  (v>-^')  } . 

Collecting  coefficients  as  indicated  in  ( 8.5)  to  (8.7),  we  get 


P - 1.  e‘"*'  { -SiHi  0^,  - e-‘*'  G , I + 

nm  J a m+l  m-l 


.iGZ  I f ( (G^j  - fG^,)  - (G^j  - fG'_j)  + 

(8.10) 

+ (nrt-l)  G . , + (n-l)  e i 1 ■*■ 

nrri  o— i. 

* S ‘ <=0.1  - <“-»  <=.-1  > * 


* < •**  <=Ll  + <=tl>  * > 
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Tte 


Inv^ 


iwvff'  I 1 + 


I '[ 

- ^ I Cl  * '■**'  Cl ' ] “„  + 


-1 


(8.11) 


+ 7^  I T ( (nri-1)  G 

iue  ‘ f 


mfl  " 


+ gL,  + 1 1 

mfl  m-1  n 


R - IT  e 
mn 


Im^ 


. fl 


M' 


(e*^  G + G , ) + 


iiH-1 


m-1' 


-1 


+ ni  I £(.‘*'  (Cl  - *'Ci>  * *■'*’  <“-i  - ''Ci>  > * 


nri-1 


(8.12) 


+ (nrt-l)  e^***  G^j^  - (m-1)  e ^'*'  1 + 


-1<|.' 


r ^ / /—  i 1 \ ^ G y 


-S-.t  ' 


nri-l  ■ ■ B-l' 

+ Cl  - ' Cl  ' ^ "n 


-l4>' 


V ■ ” I < [ “a*!  - ^ 


m-1' 


-1 


if,  c:,!  - cLl  1 ] 


(8.13) 


e‘*'G 


. (.-1)  e-‘*'  . 


+ g£.t  - G^  . 1 (fu  )'  } dz 

n+1  ®“-i- 


(8.14) 


T - - Ztr  m e 
mn 


{ inZf  G - G*  } u>  dz 
m UC.Ci  m n 


-1 
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(8.15) 


U - 2ti  e^”^'  ^ G^]  w + 

mn  J*-  ncxm-'n 

-1 


+ -^  G*  (fo)  ) • } dz 
litiC  m a 


Let  us  next  eliminate  the  Tems  which  appear  in  S 


mn 


and  U above.  Setting 
mn 


(8.14) 


J - f^7  (7  G ) (fu)  ) ’ dz 
mn  j £ £ m n 

-1 


we  have,  upon  integration  by  parts, 


1 1 


(8.15) 


J - 7 G f 0)  - (7  G„) ' fw  dz 

mn  f m n f m n 


-1  -1 


Under  the  assunption  node  on  S in  Sec.  1,  K is  bounded  on  S,  and  therefore 
it  follows,  upon  replacing  by  K in  (8.15),  that  the  first  term  on 
the  right  hand  side  of  (8.15)  vanishes,  provided  that  m.#  0 (see  (A.l)  ). 
However,  inspection  of  Q^and  S^sliows  that  we  need  never  evaluate 
if  m ■ 0.  Thus 


(8.16) 


J > f ^4-  G U)  dz  - G’  u dz,  m > 1 
mn  j f m n J m n 

-1  -1 


where  G'  - d G /dz. 
m ra 


Similarly,  we  have,  for  all  m ^ 0, 


(8.17)  (fo)^)-  dz  - - (C^)'  dz 


Some  of  the  quantities  in  (S.19)  are  not  required  for  m - o, 
since  their  coefficient  in  (8.10)  to  (8.15)  is  zero.  In  sone  cases 
this  is  fortunate,  since  same  of  the  above  integrals  do  not  exist  when 
m is  taken  to  be  zero.  In  terms  of  the  above  integrals  (8.19)  we  my 
express  the  terms(8.10)  to  (8.15)  as  follows. 


t - “M  [ e*'  I*  - e-‘''  , ) 

nri-l.n  ra-l,n  ‘ 


+ inz  I (J^  - 

nrrljii 


) - (J^  _ ) 

m-l,n  m->l,n . 


(8.20) 


+ (nri-1)  + (m-1)  e"^'^  , I 

nrfl.n 

+ — I - (m-1) 

&>£  nrH.n  m-l,n 


+ 1 } ; 
iiri-l,n  m-l,n  ’ 


Q - e*"*"  { luu  [ l‘  + e-l"  I*  , | 

™ nrfl.n  m-l.n 


- nZ  [ ] 

nrH,n  m-l,n 


8.21) 


+ [ (mfl)  (iL,  _ - 1P_,  _) 
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^ ixap  [ , , _1  , -i-f  ,1  1 

R ■ IT  e I iu»M  [ e I . , „ + e I , J 

mn  njTl,n  m— i,n 


(8.22) 


* W I * e-‘'  Wtl.a  ' '‘i-l.„> 


+ («H)  e*'  - («-»  e'"  i;.!  ,,  I 

- ^ I I^l.n  - <--l)  •’ll.n 


-iv>  ,2 


j.  -VP  1 . 

]iri-l>n  o-i»n 


-i^  ,2 


±iap  f , 1 

- TT  e { u)ul  e - e J 


(8.23) 


- ^ I (“^-1)  ^^L,n  - ""nrfl.n> 


+ (n-1)  e (I^-l.n  " ^Vl.n^ 


JP_.,  - + e'^'^  JP_,  „ 1 > ; 


Ert-l,n 


in-l,n 


(8.24) 


T = - 2tt  ine^°^  t inZ  ^ > ; 

^ ' mn  oje  mn 


(8.25) 


. 2.  ( i^U  iL  * > 


using  the  one-term  formula  (8.3),  provided  that  r'  is  not  unduly  close 
to  S . We  shall  assume  this  to  be  the  case.  VJe  shall,  however, 
reouire  G , G^,  G*,  dG  /dz,  dG^/dz  and  dG^/dz  in  order  to  evaluate  the 
integrals  (8.19).  Let  us  now  describe  the  evaluation  of  these  quantities. 

To  this  end,  let  us  set 


S { iz-z'r  + - 2f  p cos  0 }'*  , G*  5 S_^ 


m 

dG 
m 

IT 

- ^ f 

dz 

47/2 

G^  - 

'"m  1 

m 

“at/' 

3G„ 

G*  cos  md  d6 


dii  , r IK  , 

3^  - - --ir-  (-jS.  - i ) r’'  g*  cos  m9  d0 

“ 3z  .2  J R 2 


R 

d _f  if//  ^ 4.  2_x  dR  _ 1 ^ i_  j^f 


+ (—£  - i R^  ^ G*  COS  mO  d0 
R j^2^  dz 


d z 1 (if  + — •»  R*  — + ( 2.  _ i ) A 


R‘ 


+ (-^  - i )^  R*  ^ } G*  COS  me  d0 
R ^2  dz 


COS  m0  ■ 2 cos  0 cos  (m-1)  9 - cos  (m-2)  0 

then  shDW  that  given  z,z'  cind  p,  v;e  can  evaluate  the  six  quantities 
referred  to  in  the  title  of  this  section,  for  m - 0,  1,  ...  , M+1  by 
neans  of  the  following  algorithm. 


ALG0RITO1  8.1  EVALUATIOi'J  OF  G 


Evaluate  each  of  the  quantities  (19.1),  as  well  as  f - f(z)  and 

f'(z)  for  j = 1,  2,  ...  , 2N+2. 

gI,  gI,  (g  ici)',  (G^)’  ) ^ (O.O.O.O.O.O) 
in  m n m in  m 


j ^ 1 


d ■<"2c.d  d .,,in~2,  3,  ...  , Mtl 

m j m—l  m— z 

G G + G*  d 
m m j m 

G^  G^  + G*  d 

m m j J j “ 

J g’'  ■»-  G^',  + a,  G*  d 

m m ’ J j j o 

* “j  ''i  “» 

(<)'*  (<)■  + 1 + “/)  Sj  Yj  + <-j  1 CJ 

<<>'*  <8j  ♦ “/)  “j  '<i  * “j 

j : 2N+2 

<■>  «„•  <•  <■  «:>  <>' wks  <• 


9.  C0NVIR3E':CZ 


The  prxxjf  of  convergence  of  the  approxination  scheme  presented 
in  Secs.  6 to  8 is  quite  simple,  using  thm.  A.l  and  the  results  of 
Sec.  5. 

Let  us  denote  the  right  hand  side  of  (A.  36)  by  AJj^  and  for 
F e U (d,d'),  let  us  denote  .2^  by  (F) , where  (z,v>)  is 
defined  in  (5.5). 

In  view  of  Thm.  A.l,  and  Sec.  5,  we  have 

Ca)  IP„,  AJ  - Ajl„  ■>  0 for  every  Je  H (d,d'); 

(b) 

(0)8up  IP^I  < 1 + II  - <- 

Hence  according  to  t 21  pp.  469-470  J it  follows  that  the  approximation 

J • 7^  prxxiuced  by  the  algorithm  of  Sections  6 to  8 converges  to  the 
rlN 

— . 2 . 
solution  J of  Eq.  (7.2).  Moreover,  by  taking  N “ M , it  follows  that 

(9.1) 

for  some  y > 0.  Due  to  quadratures  and  matrix  soluticai  involved  in  the 
cctual  algorithm  we  actually  compute  a perturbed  solution  ho/rever, 

due  to  the  accuracy  of  the  quadrature  scheiTies  described  in  Sec.  5, 
and  since  the  resulting  Galerkin  matrix  is  not  ill-conditioned,  we 
also  have 


),  N 


i 


4 


i 

i 


(9.2) 


00 


By  combining  (9.1)  and  (9.2),  it  thus  follo'.vs  that 


(9.3)  Ij  - - 0(e"^^  ),  N 

Finally,  in  canputing  the  scattered  field  E®  » as  described  in 
Sec.  8,  we  similarly  have  by  Thm.  5.2  that 

(9.4)  I?  - . 0(e-^*■^ 

where  E denotes  the  perturbed  scattered  field  that  we  actually  conputed. 
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APPENDIX  A:  IHE  RJi4CTI0:;S  G^,  and  90^/ 3z. 

In  this  appendix  we  study  the  functions  G^,  - 3G^/3f  and 

G^  «■  3G  /3z  which  appear  in  Secs.  7 to  10.  The  results  of  this  study 
m m 

will  enable  us  to  deduce  the  following: 

i)  Each  con^xsnent  of  the  solution  j of  Eq.  (7.2)  is  iii  H(d,<i  ); 

ii)  lfia>0. 


(A.1) 


G (z.z')  - 0 (IfCz)!®) 


m 


G (z.z') 
zn 


[ G*  (z.z') 
' a 


0(1)  if  m ■ 0 

0(|f(z)r'^),  m > 0 


0(lf(z)r) 


as  z ± 1, 
z'  6 (-1.1); 


and 


(A.  2) 


(z.z')  —f 

“ f(z') 


m2 


(z.z')  - 

^ 4/f(z')(z-z’) 

G^  U.e.)  

® 4ti  f(z')(r-z') 


2'» 


as  z 
z’  e (-1.1). 


The  results  (A.l)  proved  to  be  useful  for  choosing  the  basis 
functions,  in  order  to  be  able  to  obtain  a convergent  Galerkin  method, 
while  the  results  (A. 2)  enabled  us  to  choose  the  proper  numerical 


1 

J 
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Expansion  of  the  exponential  in  (A. 6)  and  temvise 


yields 


(A.  7) 


. (-1)  V k J k K. 

1 r O r m . , o m ^ 

G — 5-  ) T^r-TT I + i o_i-i  •» 

m « 2 *• 

211  s=o 


(2s)  I 


2s+l 


where 


(A.  8) 


J®  a I { 1 - K cos^  0 f COS  (2m0)  d0 
0 

<=  _ 

K®  H I { 1 - K cos^  0}  cos  (m0)  d0 . 


The  relationship 


{ 1 - K cos^  0 } * cos  (2m0) 


(A.  9) 


{ 1 - K cos^  0}®’^  I (1-|)  cos  (2m0)  - | cos  (2ii>f2) 


shows  that 


(A.  10) 


m 


(1  - f) 


K®  - (1  - |) 

m L 


_ JS 

- J®' 

m 

4 

ntfl 

4 n* 

< 

_ <K® 

n 

" 4 

m+1 

4 m- 

whereas,  by  (A. 8) 


integration 


0-7  cos  (2m-2)0  ] 
4 


(A.  11) 


p 

i 
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The  relationships  (A.  10)  s/.ovj  that  in  order  to  eveduate  C 


m 


using  (A. 7),  we  need  only  )<now  J°  and  K°  , for  m = 0,  ± 1,  ± 2,  . . . ; 


m 


n 


we  can  then  get  the  remaining  j“  and  K for  s > 0 using  (A.IO).  To  this 

m D 

end,  integrating  temwise  in  (A. 8)  and  using  the  identity.' 


(A.  12) 


TT 

2 


cos^"  0 cos  (2m9)d0  » (^”  ) . 

DTin 


we  get 


(A.  13) 


, <y  K?“ 

^ K) 

m 2 

^ in* 


K°-  -I 

n 


IT 


If  m 


^ 0 if  m # 0 . 

In  (A.  13)  F denotes  the  hyper^eometric  function.  Ihus  we  may 


oonpute  J°  for  0 ^ k ^ . 6 by  means  of  the  formula 


(A.  14) 


,0  ~ „ 

‘'m  “ 2 ■ -2m  T (l+2ni)„  n! 

z ID*  n*u  " 


v^tereas,  if  .6  < •'  < 1 (see  [ 1,  p.559!  ] ) it  is  preferable  to  use  the 
formula 


(A.15)  I J-  - K 


1° 

m 


m 


I 

n*0 


I ) 


(i;^(n+l)  - l|/(>a+mHi)  - iiln  (1-k)  1 (I-k)" 


where 


i}/(l)  - - Y “ - .57721  56649 


(A.16)  iJ;CJs)  = - Y - 2Zn2  - - 1.9635  10026 

i|)(z+l)  - i + <j/  (z)  . 

Eq.  (A.  15)  shcjws  that 

(A.17)  J°  Hln  (1-k)  as  k -*■  T. 

m 


Now,  by  (A.  3) 


l-K  “ 


v^-4ff*  (z-z')^+(f-f*)^ 

2 “ 2 
V V 


(A.  18) 


[Itf*’^  ] (z-z*) 

4f*2 


as  z z’ 


so  that,  by  combining  (A.17)  and  (A. 18),  we  get 


(A.  19) 


j°  > in  [ ^-.T  1 , z -V  zV 


Using  induction  on  (A.  10),  it  thus  follows  that 


(A.  20) 


J®  ~0  as  z -*■  z' , s > 0 

n • 


In  view  of  (A. 7),  (A.  13),  (A.  19)  and  (A. 20)  if  therefore  follows  that 


U.21) 


G 

m 


1 fm  1 f 2f*  _ 
2tt^  (l+f*'S** 


as  z -*•  z ' , 


Si 

-I 


e, 

I 


v^ch  is  tlie  firat  of  (A.2). 

In  view  of  (A.  3),  it  follows  that 

(A.22)  K-  0(f)  as  z ± 1, 

for  all  z'  e 1-1,1  1 . Hence  by  (A. 13)  and  (A.14), 

(A. 23)  G « 0(f 1“')  as  z ■>  ± 1,  for  all  z^e  [ -1,1  1 . 

n 

f z' 

(b)  The  Rjnctions  G , G 

m ID 


Upon  defferentiating  the  expressions  (A.  3)  we  get 


(A.24) 


f _ 3v  f+f*  ^,z  = 3v  _ z-z' 
'’“3f“v  ’ 3z  V 

f ^ ^ 8ff*(ftf*) 

“ 3f  " ^2  " 

z 3k  8ff*  (z-z') 

" 3z  " ' ^3  V 


We  note  therefore,  that 


(A.  25) 


2 + V 

2 V*  + V 


X 

f 


- 0 


By  means  of  these  solutions  as  well  as  (A. 10),  we  get  the 
identities 
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(A.  26) 


(A.  27) 


^ - -1-  T f (-1)  , s-4  1 s 


4 


- iv^ko^  fl  yS+1  ^yr^S  1 \ 

2(2s+1)(28+3)  '“f  m " K V ' ^ 


- lZ5l  y { 

m 2 3 

n V 8-0 


,s.2s,  2s 


C-l)°v^°ko‘ 

(2s): 


s-1  -IN^ko^  „8  , 1 

[ (s-y  j + 2(2s+l)(2s+3) 


These  series  may  be  readily  ccsijjuted  using  the  identities 


(A.  28) 


1 TT  (3/2)„k!“ 

J ^ ” T F (^  + m,  y t m;  2m  + l;  <) 

(3/2)^ic  „ (3/2+m)  (»S+Ta) 

Z N ■«  " • k",0<k:<.6 

2 2m  , (2iiri-l)  n.  — « 

Z^ml  n-0  "^n 


=.  ^ + K (m  - h)  I 


« (4Ha)^(>Stiii) 


• I In  (l-<)  -f  2ii>  (Jstm+n)  - 2il>  (n+1)  + 1 

if  .6  < K < 1 

along  with  (A.  10).  These  identities  were  obtained  via  a procedure 
similar  to  that  vAiich  weis  used  to  get  (A.  14)  and  (A.  15). 

The  expressicwis  (A.  26)  and  (A.  28)  enable  us  to  deduce  various 
growth  properties  of  and  G*  as  z z’  and  as  z -*■  ± 1. 
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By  (A. 28)  and  (A.  3), 


(A.  29) 


(z-z')^ 


as  z -*•  z'  e (-1,1) 


and 


(A.  30) 


J ^ - 0(f*°)  } as  z ± 1 

m 


Confining  these  results  with  the  asynrototic  identities 


(A.  31) 


•c  2f* 


~ j , z ± 1, 


and 


(A.  32) 


we  get 


(A.  33) 


K 

K 


(z-z') 


z -►  z 


t 


z 

^ - 0(f)  , z ± 1 
K 


1 


z -►  z 


I 


m 


1 --  * 
4^2  ^*<*-*’)  ’ 


z' 
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Ihe  relation  (A.IO),  (A. 13),  (A. 26),  (A. 28)  and  (A. 31)  yield 


(A.  34) 


0(1)  If  D-O 
_ 0(f“"^)  if  m 


as  z ± 1, 


while  the  relations  (A.IO),  (A.13),  (A. 27),  (A. 28)  and  (A.32)  yield 


(A.  35)  - 0(f“)  as  z ± 1. 

m 

(c)  Analyticity  of  K 
Ihe  above  results  show  that 

i)  G is  bounded  as  a function  of  z on  I -1,1  1 , except  at 
n 

z ■ z'  vAiere  it  becomes  unbounded  according  to  (A. 21); 

ii)  G^  (z,z')  is  an  analytic  function  of  z e 0^,  except  at 
z - z' , Wliere  it  has  a singularity  of  the  form  (A. 21). 

In  the  following  theorem  H (d,d')  is  defined  as  in  Sec.  5. 
Theorem  A.I;  Let  K be  the  solution  of  Eq.  (4.6),  and  let  J be  defined 
by  (7.1).  Then  each  conponent  of  J is  in  H (d,d'). 

Proof ; Each  conponent  of  the  incident  field  7*  (see  Eq.  (7.3))  is  in 
H (d,d')  where  d'  > 0 is  arbitrary.  In  view  of  Eq.  (7.2),  we  need 
only  show  that  if  - (Jj^^.,  J^^)  denotes  a pair  of  functions  in 
H (d,d'),  each  of  which  is  bounded  cn  each  of  the  sets 

Sj,  - { (z,w);  z e ]a)|  - 1 } 

S2  - { (z.w)  : -1  < * < 1.  J.  < l“l  < <*’J 


Where  d*  > 1 is  arbitrary. 
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then  each  of  the  components  of 


"2  - ''2V»> 


CA.36) 


- { V 


I f ^ (nxJl)  VG 


is  in  U(d,d'). 


By  our  assun^tion  on  f in  Sec.  3»  the  solution  K ■ 

Eq.  (4.2)  is  bounded  on  S - [ -1,1  1 x [ 0 ,2Tr  J . Hence  for  z e ( -1,1  I , 
each  conponent  of  Jj^hcis  the  forro 


(A.  37) 


F (z,e^)  • I a„  (a) 


where  a /ve  H (H^);  substituting  this  fom  of  an  expression  into  (A.  36) 
m d , , 

jt 

for  J,  and  J,  and  noting  that  a /v  - 0(e  ')  V z ^ I “1»1  1 and  for 

It  1^  m 

all  d'  > 0,  we  deduce,  by  inspection  of  (7.27)  - (7.30)  and  (A.l)  and 


(A.2)  that  Ip”!.  |q"”|.  |li"|  and  |^”|  ana  0 for  all  d'  > 0 


and  for  all  z'  e [-1,1  1 . That  is,  F (z,uj)  6 H (A^,)  as  a function  of 


u),  for  all  z e I -1,1  1 . 


Hence  in  order  to  coiiplete  the  proof,  we  need  cxily  shew  that 


if  0 /v  € H (fl.),  then  each  of  the  ri^t-hand  sides  of  (7.27)  to  (7.30) 
n d 


is  in  H ((2^). 


In  view  of  the  results  of  Appendix  A.  (a)  and  A.  (b) , the 


coefficients  of  0 /v  in  the  integrals  (7.27)  to  (7.30)  are  of  the 
n 


following  three  types: 
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aCZfZ')  > a(z,z')  log  |z-z'|  , aCz.z' )/(z-z' ) , 

where  a(z»z')  is  a bounded  function  in  H (fl^)  x H (^j)*  Hence,  we 
need  only  show  that  given  g each  of  the  functions  g^,  $2  and 

gj  are  analytic  in  where 

1 

gj^(z')  - I a(2,z')  g(z)  dz; 

-1 

1 

(A.38)  g2(*')  ■ I a(z,z')  log  |z-z'|  g(z)  dz; 

-1 

1 

g3(2')  - P.V.  I g(z)  dz 

“1 

It  is  obvious  that  g,  ^ H (fl .) . 

1 a 

Next  if  z'  e {Im  z'  > 0}  then 

U 

1 

(A.39)  g*(z’)  - I flff^gCz)  dz 

-1 

is  analytic  in  this  region,  and  indeed,  by  altering  the  path  of  inte- 
gration in  (A.39)  to  the  laver  boundary  of  we  see  that  g*  is  in 
fact  analytic  in  If  wa  now  return  the  p)ath  of  integration  to  the 
intervcd  (-1,1)  and  let  Im  z'  ■*  0,  we  find  that  for  z’ ^ (-1,1), 

(A.40)  g^Cz')  “ a(z',z')  g(z')  + iri  g^  (z') 

this  expression  shows  that  since  both  g^(z')  and  a(z',z')  have  an 
analytic  extension  into  so  does  g^.  Hence  g^  is  analytic  in 


64 


Finally,  writing  82  iJ^  -he  form  of  a convergent  sum 
1 

(A.41)  82(2')  - I I U-z'l  dz 

-1  ^ 

where  the  functions  and  are  in  H(n^) , we  need  only  show  that  g* 
is  analytic  in  where 

1 

CA.42)  “ I l^-z'l  dz, 

-1 

Upon  differentiating  this  expression  carefully,  we  see  that 

\(z) 

(A.43)  g*'Cz’)  - P.V.  J dz. 

-1 

By  our  argument  involving  above,  it  follows  that  g*'  is  analytic 
i.e.,  g*  and  hence  82  is  analytic  in  This  conpletes  the  proof  of 
Theorem  A.l. 

By  assunption  for  the  case  of  finite  conductivity  of  the  body 
B,  the  function  f is  such  that  the  surface  S satisfies  Liapunov 
conditions,  in  which  case  the  surface  current  K is  bounded  on  S.  However, 
one  notes  that  the  integrals  (7.27)  - (7.30)  converge  so  long  as 
f leyvl  - o (1)  as  z ± 1,  i.e.,  so  long  as  fK  - o (1)  as  z ± 1. 

Thus  our  method  gives  answers  even  if  S is  cone-shaped  at  one  or  both 
ends,  although  the  Liapunov  conditions  are  then  violated,  and  our  results 
may  have  no  physical  significance,  except  in  the  infinite  conducting  case. 
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